Statistics of Weighted Tree-Like Networks 
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We study the statistics of growing networks with a tree topology in which each link carries a 
weight (kikj) e , where ki and kj are the node degrees at the endpoints of link ij. Network growth 
is governed by preferential attachment in which a newly-added node attaches to a node of degree k 
with rate A^ — k + A. For general values of 9 and A, we compute the total weight of a network as a 
function of the number of nodes N and the distribution of link weights. Generically, the total weight 
grows as N for A > 8 — 1, and super-linearly otherwise. The link weight distribution is predicted 
to have a power law form that is modified by a logarithmic correction for the case A = 0. We also 
determine the node strength, defined as the sum of the weights of the links that attach to the node, 
as function of k. Using known results for degree correlations, we deduce the scaling of the node 
strength on k and N. 
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I. INTRODUCTION 

The recent interest in networks stems from the discov- 
ery that the node degree distribution can have a power 
law form. Here the node degree is defined as the number 
of links that are connected to this node. It is now well 
documented that a wide variety of natural and man-made 
networks have power-law degree distributions |]|. 

In most previous studies, it has been implicitly as- 
sumed that the quality of each link is identical. However, 
there are many examples of networks in which the weight 
of each link can be distinct. For example, in scientific 
collaborations, co-authors can have a variable number of 
joint publications. Thus in the associated co-authorship 
network, a link between two individuals represents the 
collaboration strength, such as the number of jointly- 
authored publications |2|, y]. 

Transportation networks represent another situation 
where the weight of each link is variable. In the specific 
example of airline route networks, links between two air- 
ports represent the passenger capacity on this route 0. 
It has been found that the weight of the link (defined as 
the number of available seats on flights between the two 
airports connected by the link) is well-approximated by 
the products of the endpoint degrees raised to a power |5| . 
A related line of research focuses on general properties of 
flows in heterogeneous networks 0, Q • 

A third example where link weights are connected to 
network topology occurs in the metabolic network of the 
bacterium E. coli. Here a link between two metabolic 
substrates represents an enzymatic reaction, and the 
weight of a link (defined as the metabolic flux rate) dis- 



plays the same functional dependence on the endpoint 
degrees as that of the world airline network p| . 

There has also been much very recent work on var- 
ious aspects of weighted networks including character- 
izin g p aths on these networks structural properties 
0P, and the application to weighted networks to di- 
verse physical systems, such as earthquakes ^2 an d to 
synchronization phenomena |13| . 

In this work, we characterize basic properties of grow- 
ing networks in which the links possess variable weights 
|l4j . We consider the generic situation in which network 
growth is governed only by the node degrees, so that 
the link weights are passive subsidiary variables. (The 
comlementary situation in which the link weights deter- 
mine network growth has also been analyzed, see, e.g., 



Refs. 8, 9, 10 



growth 

nam 



). We define the weight Wij of a 



link between nodes i and j in the network to be the prod- 
uct of endpoint node degrees; that is, — kikj. More 
generally, we study the general case where the link weight 
has the following dependence on the endpoint degrees: 



(kikj) 6 



(1) 
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We show that this weighted network has a surprising de- 
pendence of the total weight on the system size and an 
unusual scaling for the distribution of link weights. Many 
of the new results that we obtain follow naturally from 
previous work on the degree correlations in growing net- 
works [l6|. 

For the network growth mechanism, we employ prefer- 
ential attachment, in which a new node links to previous 
nodes with an attachment rate Ah t hat depends only on 
the degree of the "target" node [Fif. For simplicity, we 
consider the situation where the newly-introduced node 
connects to a single target. Because of this topologi- 
cal restriction, the resulting network is a tree, a feature 
that turns out to simplify our analysis. To complete the 
description of the system, we need to specify the attach- 
ment rate A^. We will consider the Simon model [l8| . 
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where = k + A, with A > — 1 but otherwise arbitrary. 
For this growth mechanism the degree distribution has a 
power law tail oc k~ 3 ~ x . 

In the next section we compute the total network 
weight W(N), defined as the sum of the link weights 
Wij , as a function of the total number of nodes N. The 
distribution of link weights is investigated in Sec. Ill, and 
the node strength, defined as the sum of the weights of all 
the links that attach to the node, is studied in Sec. IV. 
We first focus on the case of the strictly linear attach- 
ment rate where Ak = k and then outline corresponding 
results for the more general case of Ak = k + A with 
A > — 1 but otherwise arbitrary. Numerical simulation 
results and comparison to our analytic predictions are 
given in Sec. V. 



II. TOTAL NETWORK WEIGHT 

We denote by Nk,i(N) the number of nodes of degree 
k that attach to an ancestor node of degree I. This de- 
gree correlation function is well-defined only for growth 
processes in which each node attaches to exactly one pre- 
vious node, so that the ancestor of each node is unique. 
Notice that the correlation function Nk,i is not symmet- 
ric with respect to the end nodes. For graphs that are 
built sequentially by preferential attachment, we are in- 
terested in the total weight of the network, defined by 



W(N) = J2klN k ,i(N). 



(2) 
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That is, we first consider the case of 9 = 1 in Eq. JIJ. 

For growth processes in which the attachment rate Ak 
does not grow faster than linearly with fc, the correlation 
function Nk,i(N) grows linearly with the total number of 
nodes N (see [l6| , and also [jjj for a brief review) . That 
is, Nk,i(N) — Nrikj for large N. The explicit form of the 
distribution nu for the simplest case of the Barabasi- 
Albert (BA) model (A = 0) is 



rihi 



4(1-1) 



k(k + 1)0 + 00 + I + 1)0 + 1 + 2) 
12(1-1) 



' k(k + l- l)(fe + 00 + 1 + !)0 + 1 + 2 ) ' 

Since W(N) — N^klnk.i, one might naively antic- 
ipate that the total weight will scale linearly with N . 
This expectation is erroneous because the infinite sum 
^2 kl nk,i formally diverges, leading to an additional N 
dependence in W(N). 

To compute ^ kl nu,i we write m = k + I and recast 
the sum as 



kln k ,i 

k.l 



= E 

m>2 

+ E 



m(m + l)(m + 2) m + 1 — I 



m— 1 

E 



i(i-i) 
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_ 2 (m - l)m(m + l)(m + 2) ^ 



For large m, the internal sums behave as 

m ~ 1 1(1 - 1) 

E— '— = m 2 h\m + 0(m 2 ) 
m+l-l 

1=2 
m — 1 



j2 1(1-1) = \m 3 + 0(r 



1=2 

Thus the asymptotically dominant part of the sum is 

E^-E^F + E^- 1 ) w 

k,l m>2 rn>2 

To evaluate these sums explicitly, we need to impose a 
cutoff in the upper limit. For the case A = 0, we there- 
fore use the fact that the maximal degree k max scales as 

n 1 ' 2 mill. This result is easy to derive by utilizing 

Nk ~ 4Nk~ 3 and estimating the maximal degree from 
the extremal criterion that in a network of N nodes, 
there will be a single node whose degree in the range 
Omax, 00 ). This criterion corresponds to the condition 
J2k>k = 1. Thus we cut off the sums in Q at 

m m ax ~ (27V) 1 / 2 to give 



rikj 



i(lnA) 2 



0(ln iV). 



(5) 



Since the leading sum in 10} diverges logarithmically, the 
mere knowledge of the scale of the cutoff is sufficient to 
determine the exact asymptotics. With this result for the 
internal sums, we thereby find, from @, 



W(N) = - A(lnA) 2 +0(NlnN). 



(6) 



A similar computation for the Simon model (where 
Ak = k + A) with arbitrary A is difficult to perform ex- 
actly because the expression for rik,i is an infinite series 
[ill rather than a rational function as in Eq. J3|). How- 
ever, the dominant contribution to the sum , kl nk,i is 
concentrated in the region 1 < < /, where n.^; admits 
the simple asymptotic form pl| 



k~ 2 l 



-2-X 



Therefore 



E^M-E'-^E^-E'- 



l-A 



(7) 



In I. (8) 



k,l 



For A > the sum converges, while for A < the sum 
formally diverges and we must again determine an upper 
cutoff. Applying the previously-mentioned extremal ar- 
gument for the maximal node degree for general values of 
A, we now find fc max oc A 1 /( 2+A ) H3. Thus we conclude 
that the last sum in Eq. © scales as 7V~ A /( 2 + A ) In N for 
— 1 < A < 0. In summary, 



W(N) oc < 



N A > 0; 

N (In A) 2 A = 0; 

N 2/(2+\) lnN A<a 



(9) 
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Consider now the situation where link weight depends 
on its endpoint degrees according to Eq. (J). We treat 
only the case 9 < 1 because this regime corresponds to 
real networks Using the asymptotic behavior for 

rik,i given in Eq. (|7|), we then have 



2-A+fl 



k.l 



which converges for A > 9 — 1 and diverges otherwise. 
Thus in the range 9 strictly less than 1, the total weight 
of the network scales according to 



W(N) ~ < 



N 

N\nN 

N {i+e)/{2+\) 



X > 6- 1; 
X = 6-l; 
X<6-1. 



(10) 



III. LINK WEIGHT DISTRIBUTION 

In this section, we consider exclusively the case of 9 = 
1. We define the link weight distribution H(w) as the 
fraction of links whose weight equals w. That is, 



nk,, 



(11) 



kl—w 



This natural definition leads to erratic behavior in the 
distribution, as the arithmetic nature of w plays a signif- 
icant role. If w is prime, then the only possible degrees of 
the end nodes are 1 and w, while if w has many divisors, 
there are more possibilities for the degrees of the end 
nodes. Thus in addition to an expected global smooth 
dependence of IL(w) on w, there will be superimposed 
fluctuations related to the divisibility of w. These fluc- 
tuations will disappear if we average H(w) over a range 
much less than w and much larger than the average dis- 
tance between successive primes (which grows as hiw 
pnjl- This suggests that we also consider a smoothed 
distribution V{w) defined by 



E 

■w<.kl<w-\-w a 



n>k,i, 



(12) 



with a in the range < a < 1. For sufficiently large 
w this smoothed distribution should be monotonic and 
independent of a. 

To determine the behavior of the link weight distribu- 
tion for id > 1, we again start with the case of the BA 
model where A = 0. Using J2J and writing I — w/k we 
obtain, for large w, 



IL(w) 



E 

k\w 



I2wk 2 



(k 2 + w ) 3 {k 2 + w y 



(13) 



The notation k\w means that this sum runs over all values 
of k that are divisors of w. For large w, the asymptotic 
behavior of the summand scales as w~ 2 times the average 



number of divisors of w. Using the celebrated Dirichlet 
formula (see e.g., HJ), the total number of divisors of all 
integers from 1 to w is 



d[x) — white + (27 — l)w + 



where d(x) is the number of divisors of the integer x and 
7 = 0.577215 is the Euler-Mascheroni constant. Thus 
the average number of divisors of w grows logarithmi- 
cally with w; more precisely, by differentiating the above 
formula with respect to w, we have (d(w)) — > law + 27. 
Thus the tail of the weight distribution is 



V(w) 



hi 1 



as w ^> 1. 



(14) 



Consider now the general case of arbitrary A (but with 
the necessary constraint A > — 1 to ensure a finite average 
degree). In this case, we have to combine two separate 
contributions to determine the link weight distribution. 
First, using kl = w, we re-write Eq. as 



-2-A 



(15) 



This asymptotic behavior holds for k < I, so it can be 
used to estimate the part of the sum in Eq. (|11|) that runs 
over the divisors k\w such that k < y/w. The correspond- 
ing contribution to the weight distribution is therefore 



V<(w) - w- 2 - x ax(Vw~) 



(16) 



where o\ (ri) is the generalized divisor function of n de- 
fined as 



/In 



(17) 



where the sum runs over all divisors / of the integer n. 
This generalized divisor function exhibits, on average, 
a power-law growth in n for A > 0, cr\(n) ~ n x , and 
approaches a constant for A < [22]. As a result, the 
asymptotic behavior of V<{w) undergoes a qualitative 
change when A passes through zero. Using the result 
for the generalized divisor function in Eq. (|16fl . we find, 
for the contribution to the link weight distribution from 
divisors of w that are smaller than \/w: 



r<(w) 



-2-X/2 A > 0j 



A < 0. 



-2-A 



(18) 



A second contribution to the weight distribution arises 
from the range k > I. (For the following arguments we 
use the smoothed weight distribution of Eq. (|12|) . Corre- 
spondingly, the generalized divisor function a\ (n) should 
also be averaged in the same way as that for the weight 
distribution in Eq. (|12|) .) In this case, the asymptotic 
behavior of the distribution nk,i is |16| 



n k ,i ~ I 



6+3A 



1 1' 



-5-2A 



(19) 
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We use this asymptotics to estimate the part of the sum 
in Eq. that runs over the divisors k\w such that 

k > y/w, or cquivalently over divisors l\w such that 
I < \fw. The corresponding contribution to the weight 
distribution is 



(20) 



Since 6 + 3A > (in fact, the stronger inequality 6 + 
3A > 3 holds), the generalized divisor function behaves 

as cr 6+3A (Vw) - w 3+3A/2^ Thug 



-2-X/2 



(21) 



We therefore find that the large-w tail of the complete 
weight distribution V(w) — V<(w)+V>( , w) is dominated 
by 'P< (w) . Thus summarizing our results, the link weight 
distribution is 



V(w) 



w -2-X/2 X>Q . 

2 \nw A = 0; 
- 2 - A A<0. 



(22) 



The average weight per link (w) — J dw w Tl(w) should 
be equal W(N)/N. In computing (uu) using Eq. I)22[l. we 
find agreement with Eq. @ only for A > 0. The reason 
for the disagreement in the range — 1 < A < is unclear. 
Part of difficulty may lie in the fact that the link weight 
distribution exhibits an anomalous feature that is not 
simply a power-law behavior (see Sec. V). 



IV. NODE STRENGTH 

Many real networks display a significant correlation be- 
tween the degrees of adjacent nodes. There are examples 
where there is a bias for nodes to be connected to other 
nodes with either similar (assortative) or different (dis- 
assortative) degree j23|. This tendency to select nodes 
with a consonant property is often expressed as a Pear- 
son correlation coefficient of the network. For networks 
where the weights are correlated with the topology, we in- 
vestigate the effects of correlations by studying the node 
strength defined as Q 



Si 



(23) 



where the sum runs over all nodes j that are linked to 
node i. Let us first consider the simplest link weight func- 
tion where Wij is just the product of the endpoint node 
degrees. Then the strength s(k) of a node of degree k is 
simply k times the sum of the degrees of the neighboring 
nodes. Since there are k neighbors, s(k) is then k 2 times 
the average degree of the neighboring nodes if we ignore 
correlations between the degrees of neighboring nodes. 
Now the average node degree equals twice the number of 
links divided by the number of nodes; this simply equals 



2 for a network with a tree topology. Hence the node 
strength s(k) scales as 



s(k) ~ k 2 



(24) 



Because node degrees are actually correlated, we now 
investigate their influence on the simple prediction of 
Eq. (J23J. To include the effect of node correlations on 
the behavior of the node strength s(k), we start by split- 
ting the contributions to s(k) into two parts: one from 
incoming links and another from the outgoing link. This 
gives 

s(k) = Sin(fc) + s ou t(fc) = (k-l)k(jk)in + k(h)out- (25) 

The first term accounts for the fact that there are k — 1 
incoming links, each of which has weight kj k , where jk 
denotes the degree of a daughter node when the initial 
node has degree k. Similarly, the second term accounts 
for the single ancestor node with degree Ik- We now 
consider these two contributions to s(fc) separately. 

For the contribution due to the ancestor node, we use 
the fact that its average degree (ifc) ou t is given by 



(I 



k) out 



(26) 



E;>i n k,l 

The normalization factor is just the degree distribution 
x 4 



nk ,i = n k 



k(k + l)(k 



Using this, together with Eq. ©, we obtain (again using 
the shorthand notation m = k + I) 



(Uout = (fc + 2)^ 



1(1-1) 



i>i 



m(m + l)(m + 2) 



+3(fc + l)(fc + 2)^ 



1(1-1) 



(m — l)m(m + l)(m + 2) 



The first sum on the right-hand side is logarithmically 
divergent. Using the fact the upper cutoff is v^V and 
replacing the sum by an integral, we find that for k 3> 1 
the first sum is asymptotic to In y/N — In k. Similarly 
we find that the second sum is asymptotic to l/(3fc). 
Therefore 

Sout(fc) = k(l k ) out -> k 2 (lnViV-lnfc) . (27) 

For the contribution to s(k) from the daughter nodes, 
we need the average degree of these nodes. This is given 
by 



{jk 



Sj>1 n j,k 

The denominator is found from 

4(fc-l) 



(28) 



^ n jt k = (k- l)n k = 



fc(fc + l)(fc + 2) 
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This equality is somewhat subtle. A node of degree k has 
k — 1 incoming daughter nodes. If we sum the two-point 
correlations functions for all these daughter nodes, we 
overcount the number of nodes of degree k by exactly the 
factor k— 1. The numerator can be computed straightfor- 
wardly but tediously by following similar steps to those 
preceding Eq. I|27|) to evaluate the sums. The result- 
ing asymptotic behavior has the simple form 4fc~ 2 lnfc. 
Thus we have 

Si„(fc) = (k - l)k O fe )in -> k 2 lnk (29) 

Consequently, for k ^> 1 the node strength scales 

s(jfc) -fc 2 lniV. (30) 

A note of warning is in order. The above compu- 
tation of the contribution of the outgoing link to the 
node strength is exact — both s ou t(k) = k {Ik} out and 
Eq. I|26[> are manifestly correct, and therefore the knowl- 
edge of njy suffices to determine s ou t(fc). On the other 
hand, an exact calculation of the contribution of incom- 
ing links requires knowing the many-body correlation 
function n(ji, . . . ,j k -i)- This quantity gives the fraction 
of nodes of degree k whose daughter nodes have degrees 
ji, . . . ,jk—i- This exact contribution has the form 

S in (fc) = k (ji + . . . + jk-l)m , 

Eji,...,^-! (ii + • • • + 3k-i)n(ji, . . .,jk-i) 

E 3 - 1 ,...j»_ i n(ji»-" J ifc-i) 

Unfortunately, we do not know n(ji, . . . , jk-i), so we 
cannot compute S[ n (k) exactly. Equations and (J2BJ, 
which involve only the known degree correlation n^z, are 
based on the assumption that the degrees of the daugh- 
ter nodes are uncorrelated. Although the correlation be- 
tween the degrees of daughter nodes is certainly smaller 
than the degree correlation between the daughter and 
mother nodes, it is not clear that we can ignore these 
correlations. 

Consider now the general linear attachment kernel 
Ak — k + A with A > — 1 and choose the weight func- 
tion Wij = (kikj) e with 8 < 1. We now compute the 
contribution of the outgoing link to the node strength. 
We have 

k e x 

sout(fc) - — }J n kt i (31) 
nk 

i>i 

In performing this sum, we use the asymptotics ~ 
/c~ 3 ~ A for k 3> 1 and the results of [16j for the asymp- 
totics of nk,i, 

n k ,i ~ k~ b - 2X l 1+x for l<Z<fe 

nk,i ~ k~ 2 l~ 2 ~ x for l<fc<i. 

We estimate the formally diverging sums by using the 
cutoff 7Vi/(2+A)^ We thus find 

(k 26 \>9-l; 
flout (A) ~ { k 2+2X (\n 7VV(2+A) _ lnfc) \ = 6-l; 

[ k i+\+e N (o-i-\)/&+x) _ k 2e x<6-l. 



Finally, we specialize this prediction to the simplest 
weight function my = as we report simulation re- 

sults only in this particular case in the next section: 

(k 2 A > 0; 

flout (*0 ~ < fc 2 (lnViV-hifc) A = 0; (32) 

[fc2+A iV -A/(2+A) _ k 2 X<Q 

V. SIMULATION RESULTS 

We now present numerical results for the total weight 
of a network and the underlying weight and strength dis- 
tributions. To generate the network, we use the redirec- 
tion method, as discussed in Ref. [ifj. In this approach, 
we attach a new node to a randomly-selected target node 
with probability 1 — r and attach to the direct ancestor 
of this target with probability r. This method is both 
extremely simple, because the target node is randomly 
selected, and efficient, since the time to build a network 
of N nodes scales linearly with N. 

As discussed in Ref. , this redirection rule is equiv- 
alent to a network in which the attachment rate Ak to 
a node of degree k equals k + A, with A = - — 2. The 
attachment rate k + A then leads to a degree distribution 
nk ~ with v = 3 + A 16]. Thus by implementing 
redirection, we can generate scale-free networks with de- 
gree distribution exponent anywhere in the range (2, oo). 
As a parenthetical technical point, each node in the initial 
network must have a unique ancestor to ensure that redi- 
rection weights each network realization correctly po| . In 
our simulations, we chose the starting network to be a tri- 
angle in which each node points cyclically to its nearest 
neighbor. This detail about initial condition does not 
affect the asymptotic form of the degree distribution. 

Figure ^ shows the average network weight as a func- 
tion of the total number of nodes N. These results are 
based on 100 realizations of the network that each con- 
tain up to A^ max = (1.5) 40 w 1.1 x 10 7 nodes. On a double 
logarithmic scale, the data for W(N) versus N is quite 
linear for A > 0, but there is a very small systematic cur- 
vature in the data for A < 0. As suggested by Eq. ©, we 
therefore divide the data for W(N) by (In AQ 2 for A = 0, 
or by In N for A < 0. 

The resulting data shown in Fig. |2 are now visually 
more linear on a double logarithmic scale. Power-law fits 
to large-A subranges of the data give effective exponents 
that are only weakly dependent on the subrange. For 
the case A = —2/3, the effective exponent changes from 
1.451 to 1.487 as the lower limit of the fit range iV m i n is in- 
creased from 4 to 4 x 10 5 while A max remains fixed. This 
behavior suggests that the total network weight scales as 
AT 3 / 2 In A, as predicted by Eq. ©. Similarly, for A = 0, 
the effective exponent changes from 0.938 to 0.969 when 
the above fitting protocol is applied. Again, the data is 
consistent with W(N) ~ N (In A) 2 for A = 0. For A > 0, 
the variation in the effective exponent is smaller still and 
it is evident that W(N) varies linearly with N . 
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FIG. 1: Average network weight W(N) as a function of iV for 
the representative cases A = 100, 1, 0, and —2/3, (o, +, A, V), 
respectively. 
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FIG. 2: Same as in Fig.0 except that the data for A = and 
A < has been divided by (IniV) 2 and by IniV, respectively. 
The solid lines are best fits to the data for N > 4 x 10 5 . 
The slopes of these lines are 0.990, 0.988, 0.969, and 1.489 for 
A = 100, 1,0, and —2/3, respectively. 



In Fig. [2 we show our data for the distribution of link 
weights for the representative cases of A = 0, 1, and —0.5. 
The data shown are based on 10 4 network realizations 
with TV = 10 6 for A = and 1, while N is limited to 
25,000 for the A = —1/2 due to the broadness of the 
link weight distribution. The data shows considerable 
small-scale variation because of the previously-mentioned 
fluctuations in the number of divisors of the integers; 
thus we also examined the locally averaged distribution 
given by Eq. (|12|) . While this construction smooths the 
distribution, exponent estimates based on the smoothed 
distribution are nearly identical to those found for the 
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FIG. 3: Distribution of link weights U(w) versus w for A 
(top), A = 1 (middle) and A = —0.5 (bottom). Every 
tenth point is plotted and the line is a 10-point average. The 
quantity II(ui)/lnui appears in the top plot. 



raw distribution (H(w) in Eq. Ijllfl ) and thus we quote 
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FIG. 4: The node strength s(k) versus k for a network of 
N = 10 7 nodes with A = 0, and the contributions to this 
strength from incoming and outgoing links, as given by the 
two terms in Eq. 1251 (□,«, A respectively). Inset: s(fc)i n 
divided by k 2 In k. 



results based on the analysis of the latter. 

For A = 0, a power-law fit to the data for H(w)/\nw 
gives n(ti;)/ lnu> ~ w~ x , with x ~ 2.3, compared to the 
theoretical prediction from (|22|l of x = 2. For A = 1, we 
find H(w) ~ w~ x , with x ~ 2.23, while the theoretical 
prediction from (|22|l is \ — 2-5. Finally, for A = —0.5, 
a simple power law fit to the raw distribution in Fig. |3| 
is clearly an oversimplification. However, a power-law fit 
to all the data gives x = 1-39, close to the theoretical 
prediction of x — 1-5. Overall, the agreement between 
the predictions of Eq. (1221) and simulation results is sur- 
prisingly good, given the vagaries of these distributions. 

We also numerically calculate the node strength as 
function of the degree k. In Fig. 01 we show the results 
for A = for a network of N = 10 7 nodes. As a check of 
Eq. (J37J, the contribution to the node strength from the 
outgoing links was fit with form s ou t{k) ~ fc^(lnfc max — 
lnfc). Using a non-linear curve fitting package, we ob- 
tain ft = 1.9 and fe roax = 13,000. This value of fc max 
is very close to what we observed numerically (10,900) 
for the maximum degree over all 1000 realizations of the 
network. Overall, the data for s ou t(fc) agrees reasonably 
well with the theoretical prediction from Eq. H27|) with 
= 2. 

In Fig. we show the node strength for two other rep- 
resentative values of A. For the case of A = 1, both 
s(k) and s out (fc) grow as fc 2 , as predicted by Eqs. l(2"Y|) 
and l)3()[l. For the case of A = —2/3, s ut(fc) grows as 
as predicted by (|2*7|) . while Si n (k) grows slightly 
faster than k 2 , which is consistent with the k 2 In k growth 
given in Eq. (|29ll . Notice also that the point where Si n (k) 
and Sout(^) intersect moves to larger k as A decreases. 
This crossover point also coincides with /c max and the 
behavior of fc max on A is consistent with the prediction 
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FIG. 5: The node strength s(k) and its contribution from in- 
coming and outgoing links (□, •, A respectively) for the cases 
of A = 1 (upper) and -2/3 (lower), both with N = 10 7 . 



While most of our discussion has considered the case 
where the exponent 9 = 1, there are examples of real 
networks, such as the E coli metabolic networks Q and 
the airline route network [5j, where the correlations be- 
tween node degree and link weight seem to follow Wij = 
{kikj) 1 / 2 , that is 9 — 0.5. To test the predictive power 
of our analytical approach also for the case of 9 < 1, we 
investigate the scaling of node strength as a function of 
node degree for 9 = 1/2 for the two specific examples of 
A = and A = —2/3 (Fig. EJ). Our analytical results for 
the outgoing node strength are s ou t{k) ~ k (A = 0) and 
s out(&) ~ k 5 / 6 (A = —2/3), in good agreement with the 
best fits of s out (fe) ~ k 9 and s out (fc) ~ k°- 8 , respectively. 

As a further test of our theoretical predictions, we show 
data for the node strength in which we attempt finite-size 
scaling (Fig.EJ. In panel (a), the data for A = is consis- 
tent with the asymptotic behavior of s ou t (k)/ In N ~ k 2 
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FIG. 6: The node strength for 6 — 0.5. The left panel shows 
the case A = and the right A = —2/3, both for networks of 
size N = 10 6 . The node strength s(k) and its contribution 
from incoming and outgoing links are denoted by □, • and A 
respectively. 



given in Eq. (|27(l . The departure from data collapse arises 
from the correction term in this equation. In panel (b), 
the data for A = 2/3 is now consistent with the form 
SoutOV^ 172 ~ fc 4/3 - k 2 jN x l 2 given in Eq. fgj. The 
departure from data collapse is due to the influence of the 
correction term in s ou t(k). Panel (c) shows the expected 
scaling behavior of S[ n (k) ~ k 2 \nk given in Eq. (|29|l . 
Curiously, this same behavior also occurs for the case of 
A = -2/3. 



VI. SUMMARY 

We studied the statistics of growing networks in which 
each link ij in the network has an associated weight 
(kikj) 6 , where fcj and kj are the node degrees at the end- 
point of the link ij. We also characterized a node by its 
strength, defined as the sum of the weights of the links 
that are attached to this node. The link weights and node 
strengths provide new metrics with which to characterize 
heterogeneous networks. The motivation for considering 
these properties stems from a number of real-world ex- 
amples, such as airline route networks or metabolic net- 
works, in which the flow capacity of each link is different. 

We focused on developing an analytical understand- 
ing of the link weights and node strengths as a func- 
tion of the underlying growth mechanism of the network. 
Much of our analytical results are obtained by adapting 
previously-derived results for node degree correlations to 
the link weights and node strengths. Generically, we 
found that the total network weight grows faster than 
linearly with the total number of nodes N. Strictly lin- 
ear growth of the total weight on TV occurs for networks 
whose degree distribution decays relatively quickly in k 
for large k. The distribution of link weights has a power- 
law tail that is modified by a logarithmic correction for 
A = 0. There is also a strange and, as of yet, unexplained 
anomaly in this distribution when A < 0. There are also 
small-scale fluctuations in the link weight distributions 
that are caused by the granularity in the number of ways 
that a given link weight can be factored into a product 
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FIG. 7: Finite size scaling of the node strength contribution 
from outgoing ((a) and (b)) and incoming ((c) and (d)) links 
for A = (upper panels) and A = —2/3 (lower panels). The 
network sizes are N = 10 4 , 10 5 , 10 6 and 10 7 (□, •, A, * respec- 
tively) . 



of the degrees of the endpoint nodes. 

A calculational approach in the same spirit as that used 
for the link weights was also applied to determine the 
node strength. For A > 0, the node strength scales as k 2 , 
with a logarithmic correction in the case of A = 0, while 
for A < 0, the node strength has a more complicated scal- 
ing behavior. The node strength can also be decomposed 
into contributions from the outgoing link and a contribu- 
tion from incoming links. Only the former is amenable to 
a complete analytical treatment in terms of the two-node 
correlation function, while the latter requires knowledge 
of many- node correlations. We expect that these many- 
node correlation functions can be computed withing a 
rate equation approach. 
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APPENDIX A: LOWER BOUND FOR THE 
TOTAL LINK WEIGHT 

Our calculation for the total network weight relied on 
the distribution nk,i that was obtained in Ref. [lq|. We 
also made use of the cutoff fc max that applies for the de- 
gree distribution n/- pfij . Thus the application of this 
cutoff to rik,i could in principle be questioned. Here we 
establish (for the simplest A = model) a rigorous lower 
bound for the total weight of a network using only the 
well-known properties of the degree distribution. 

Let W(N) be the average total weight of the network 
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with N nodes. The next node will link to a target node of 
degree k with probability | kn k . The degree of the target 
node then increases k + 1, so that the average weight of 
the new link is ^2(k + 1) x (krik)/2. The increase of the 
degree of the target node (k — ► fc + 1) implies that the 
weights of the k current links that end at the target node 
will also increase; we denote this corresponding average 
increase by (...)• Thus 



W(N + 1) = W(N) +J2l k ( k + + (• 



(Al) 



We shall ignore the positive increment (...) since its cal- 
culation requires knowledge of the distribution nu,i- By 
this approximation we will derive the lower bound rather 
than the true asymptotic. 

Equation (|A1|) . in conjunction with the well-known ex- 
pression for the degree distribution 



fc(fc + l)(fc + 2) 



then gives 



W(N + 1) > W(N) + Y^ 



k + 2 



(A2) 



(A3) 



For the degree distribution, the location of the cutoff 
fcmax ~ N 1 / 2 is well established. (Even the full scal- 
ing form nfc(iV) = n k F(k/y/~N) is known analytically 
pj.l Thus W(N+1)-W(N) > hxN, so that the lower 
bound for the total weight is 



W(N) > N\nN. 



(A4) 



This bound exhibits a slower growth than the prediction 
of Eq. © but it proves that the total weight of a network 
grows faster than linearly in N. 

For the general case where the attachment rate as the 
form Ak = k + A, the tail of the degree distribution is 
given by n k ~ fc~ 3 ~ A . For A > this asymptotic form for 
the degree distribution implies that the sum k(k+l)nk 
converges, leading to the lower bound W(N) > N. This 
gives the same scaling behavior as the actual asymptotic 
behavior of Eq. QJ). For A < 0, the sum ^ k{k + l)nu 
diverges and we use the cutoff fc max ~ N 1 /^^ to obtain 
J2k(k + l)n k ~ AT-V(2+A). This then 

gives the lower 

bound W(N) > N 2 ^ 2+x \ a result that also agrees with 
the actual asymptotic behavior in Eq. ©. 
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